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■ We calculate the cross section for the inclusive production of single hadrons 
with finite transverse momenta in deep-inelastic scattering at next-to-leading or- 

| der (NLO), i.e. through O(a^), in the parton model of QCD endowed with non- 

perturbative parton distribution functions (PDFs) and fragmentation functions (FFs). 
The NLO correction is found to produce a sizeable enhancement in cross section, 
i-£h I of up to one order of magnitude, bringing the theoretical prediction to good agree- 

i" 1 ■ ment with recent measurements for neutral pions and charged hadrons at DESY 

HERA. This provides a useful test for the universality and the scaling violations of 
the FFs predicted by the factorization theorem. Such comparisons can also be used 
to constrain the gluon PDF of the proton. 
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1 Introduction 



The predictive power of the parton model of quantum chromodynamics (QCD) lies in 
the factorization theorem. In deep-inelastic scattering (DIS), factorization in short- and 
long-distance parts allows us to describe the observed cross sections of inclusive hadron 
production as a convolution of the partonic cross sections with non-perturbative parton 
density functions (PDFs) and fragmentation functions (FFs) pQ. Single-hadron inclusive 
production in electron-proton DIS, 

e-(k)+p(P)^e-(k') + h(p)+X, (1) 

occurs partonically already in the absence of strong interactions, at 0(ot®), where a s is 
the strong-coupling constant, when one parton of the proton (a quark) interacts with the 
lepton current and fragments into the hadron h (naive parton model). If the virtuality 
Q 2 = —q 2 of the four-momentum transfer q = k! — k satisfies Q 2 <C m|, where mz 
denotes the Z-boson mass, then process (0) is essentially mediated by a virtual photon 
(7*), while the contribution from Z-boson exchange is negligible. In the following, this is 
the situation we are interested in. 

Since we are interested in perturbative QCD effects, we require the hadron to carry 
non-vanishing transverse momentum (p^) in the centre-of-mass (cm.) frame of the virtual 
photon and the incoming proton. At leading order (LO), the corresponding partonic 
subprocesses thus contain two partons in the final state, one of which fragments into 
the hadron, while the other one balances the transverse momentum. At next-to-leading 
order (NLO), three-parton final states contribute to the real correction, while the virtual 
correction arises from one-loop diagrams with two final-state partons. 

The investigation of single-hadron production is interesting for several reasons. First 
of all, it provides a test of perturbative QCD and of factorization. Apart from the partonic 
cross sections obtainable from perturbative QCD, the theoretical predictions essentially 
depend on universal PDFs and FFs. In particular, the FFs, which are fitted to electron- 
positron-annihilation data, may be tested with regard to their universality. Furthermore, 
the theoretical predictions allow for a direct comparison with experimental data, without 
resorting to any kind of Monte Carlo model to simulate the hadronization of the outgoing 
partons. Thus, we may expect very meaningful results. Moreover, the theoretical predic- 
tions are directly sensitive to the gluon PDF of the proton with the potential to constrain 
the latter. 

On the experimental side, precise data were collected by the HI |2|3|4j and ZEUS |5l6j 
Collaborations at the ep collider HERA at DESY. They refer to 7r° mesons in the forward 
region |3l4j . with small angles with respect to the proton remnant, and to charged hadrons 

mm- 

More than 25 years ago, the cross section of process (JTJ with finite transverse mo- 
mentum of the hadron h was calculated by Mendez [T at LO, to 0(a s ). Since QCD 
corrections are typically large and we are confronted with precise experimental data, it is 
desirable to compare these data with predictions of at least NLO accuracy, including the 
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terms of 0(a 2 ). For this purpose, a first NLO QCD prediction was computed in Ref. jH], 
neglecting the longitudinal degrees of freedom of the virtual photon. 

The theoretical description can be rendered more reliable by resumming the leading 
logarithmic contributions of the perturbation expansion. In this sense, the Balitsky-Fadin- 
Kuraev-Lipatov (BFKL) [9 and Dokshitser-Gribov-Lipatov-Altarelli-Parisi (DGLAP) [TT1] 
equations resum at LO the (a s ln(l/xs)) n and (a s ln(Q 2 / Ql)) n contributions, respectively, 
where xb is the Bjorken variable and Qq is the cut-off scale for the perturbative evolution. 
Disregarding the fact that these resummations are just approximations, they evidently 
fail in the kinematic regions of large xb values and small Q 2 values, respectively. 

In this paper, we perform a full NLO QCD calculation, also taking into account the 
longitudinal degrees of freedom of the virtual photon. We encounter ultraviolet (UV) 
and infrared (IR) singularities, which we all regularize using dimensional regularization. 
In order to overcome the difficulties in connection with the IR singularities emerging 
in different parts of the NLO correction, we employ the dipole subtraction formalism 
[TT] . In contrast to the more conventional phase-space slicing method, there is no need 
to introduce any unphysical parameter to cut the phase space into soft, collinear, and 
hard regions. Moreover, all cancellations of IR singularities occur before any numerical 
phase-space integration is performed. We thus conveniently obtain numerically stable 
predictions. 

An independent calculation was recently presented in Ref. [12 , where the matrix 
elements of the hard-scattering processes were adopted from the DISENT program pack- 
age In Ref. [12], the phase-space slicing method was applied to handle the IR 
singularities. Another related work, focusing on fracture functions of the proton, was 
published in two parts, related to incoming gluons [Hj and quarks [To] . 

This paper is organized as follows. In Section |2J we describe our analytical analysis. 
The LO result and a specific part of the real NLO correction are relegated to Appendices lAl 
and El respectively. In Section El we present our numerical results. Our conclusions are 
summarized in Sec. EJ 



2 Analytical analysis 

According to the factorization theorem P , the differential cross section for process (JTJ) is 
given as a convolution of the hard-scattering cross sections da ab with the PDFs of the 
proton and FFs D% of hadron h, as 



dx dy dz 



d 4 o h v-^ f dx f dz , . d 4 a ab . 



X 



where /Zj and /ij are the factorization scales related to the initial and final states and the 
sum runs over all tagged initial- and final-state partons, a and b, respectively. As usual, 
the dimensionless variables x, y, and z are defined as x = Q 2 / (2p a -q), y = p a -q/p a -k, and 
z = Pa-Pb/Pa-<1 with respect to the partonic four-momenta p a and pb, and their barred 
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counterparts x = xb = Q 2 / (2P ■ q), y = y = P ■ q/ P ■ k, and z = P ■ p/ P ■ q with respect 
to the hadronic four-momenta. We have Q 2 = xsyS, where S = (P + k) 2 is the square 
of the ep cm. energy. It is convenient to describe the kinematics in the cm. frame of 
the virtual photon and the incoming parton a as is done in Fig. ^ where we take the 
three-momentum of the virtual photon to point along the z axis and the three-momenta 
of the incoming and scattered electrons to lie in the x-z plane. Then, the azimuthal angle 
<p of the hadron h is enclosed between the plane spanned by the three-momenta of the 
incoming and scattered electrons and the one spanned by those of the virtual photon and 
the outgoing parton b. 




Figure 1: Cm. frame of the virtual photon and the initial-state parton a, where the 
three-momenta of the leptons are rotated into the x-z plane. 



The hard-scattering cross sections may be written as contractions of a lepton tensor 
with hadron tensors H^, as 



d 4 a ab a 2 y 

dxdydzdcf) 16n 2 Q 4 



^Ht, (3) 



where a is Sommerfeld's fine-structure constant. If the virtual photon and the initial-state 
parton are both unpolarized, then there cannot be any dependence on the azimuthal angle 
<p. Integrating over the latter, we find the decomposition 

d^ ab a 2 (y 2 - 2y + 2 b n y 2 - 6y + 6 A 

2 #t+ 2 -3^ H L ), (4) 



dxdydz 8ir \ 2yQ 2 y 3 s 
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where Hf = -g^Hfi, Hf = p%p u a Hfi } and s = (p a + q) 2 . 
The partonic subprocesses contributing at LO are 



7* + q^q + g, 
7* + q -»• g + q, 
7* + g^q + q, 



(5) 
(6) 
(7) 



where it is understood that the first of the final-state partons is the one that fragments 
into the hadron h. Here, q = qi,q l , . . . , q nf ,q nf i where n/ is the number of active quark 
flavours, which are ordered according to their masses, i.e., q\ = u, q2 = d, q-s = s, 
q4 = c, and q$ = b, and we identify q = q. There are two Feynman diagrams for each of 
the processes ©-(J7J). The matrix elements of processes ©-© are interrelated through 
crossing symmetry. Owing to charge-conjugation (C) invariance, the counterparts of 
processes ©-(J7J) with quarks and antiquarks interchanged yield equal cross sections and 
do not have to be calculated separately. However, this is not generally true for the PDFs 
and FFs. Therefore, we have to explicitly sum over all possible pairings of partons a and 
b in Eq. (J2J). For the reader's convenience, the LO expressions for the Lorentz invariants 
and Hf^ in Eq. (j3J) are listed in Appendix [XJ They are of 0(a s ) and proportional to 
e 2 , where e q denotes the electric charge of quark q in units of the positron charge. 

In order to determine the NLO correction to the cross section of processes we 
have to compute the virtual and real corrections of 0{a 2 s ) to the hadron tensors. We then 
encounter a rather involved pattern of singularities. All these singularities are regularized 
using dimensional regularization with D = 4 — 2e space-time dimensions yielding poles 
in e in the physical limit D — > 4. The integrations over the loop four-momenta in the 
virtual correction lead to UV and IR singularities, where the IR ones comprise both soft 
and collinear singularities. All UV singularities are removed through the renormalizations 
of the wave functions and the strong-coupling constant. The remaining soft and collinear 
singularities cancel partly against counterparts originating from the phase-space integra- 
tion of the real correction. The remaining collinear poles have to be factorized into the 
bare PDFs and FFs so as to render them finite. 

The virtual correction is obtained as the interference of the Born and one-loop matrix 
elements. The latter receive contributions from self-energy, triangle, and box diagrams. 
These involve two-, three-, and four-point tensor integrals, which are reduced to scalar 
integrals via tensor reduction 16J. The scalar integrals contain both UV and IR sin- 
gularities. They are computed analytically in dimensional regularization. Our analytic 
expressions for the contractions of the resulting hadron tensors with g^ v agree with the 
literature ^7j. The virtual correction is renormalized in the modified minimal-subtraction 
(MS) scheme and thus UV finite. 

The partonic subprocesses contributing to the real correction read 



l* + q^q + g + g, 
l* + q-+g + q + g, 
7* + g^q + q + g, 



(9) 
(10) 
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i* + g^g + q + q, (ii) 

7* + q^q + q + q, (12) 

7* + g->g + g + g, (13) 

7* + g - q + q + g', (14) 

7* + g-*g' + g' + g, (15) 



where q, q' = qi, q x , . . . , q nf , g with q ^ q' . As in processes ©-(J7J), the first partons in 
the final states of processes (J%|)-()15p are taken to fragment into the hadron h. The order 
in which the residual final-state partons appear is irrelevant. There are eight Feynman 
diagrams for each of the processes (J%j)-(fH?j) and four ones for each of the processes (fPfj) and 
(|15|). Crossing symmetry interrelates the matrix elements of processes flS])- (jll)) . those of 
processes (fT2 |) - (fT3|) . and those of processes (H3J)-([T5J). The cross sections of processes (JSJ)- 
(fTo]) are of (D(a 2 ). Those of processes (l%j)- (fT3j) are proportional to e; 2 while, at first sight, 
those of processes (fT^ - (fT5]l contain pieces proportional to e 2 , e q e q >, and e 2 ,. However, 
in the case of process (|T4^) . the piece proportional to e q e q r vanishes by Furry's theorem, 
as is explained below. The squared matrix elements of processes (jHj)- (jllj) involve one 
quark trace, those of processes (fT2|) and (fTHJ) contain pieces with one or two quark traces, 
and those of processes (|14|) and (J15)) involve two quark traces. Due to C invariance, the 
counterparts of processes (J%j)- (fi3j) with quarks and antiquarks interchanged yield equal 
cross sections. Notice that, in the case of process (JT^), we have to distinguish between 
the case where the tagged quarks q and q' are both particles or anti-particles and the 
case where one is a particle and the other one is an anti-particle. Processes flSJ) and (fT3|) 
each contain two identical untagged partons in the final state, so that their cross sections 
receive a statistical factor of 1/2 to avoid double counting in the phase space integration. 
We derived the matrix elements of processes dHJ (HSJ) in two steps. First, we calculated 
the ones of the corresponding processes of e + e~ annihilation via a virtual photon, which 
may also be found in Ref. [T^]. Then, we employed crossing symmetry. The squared 
matrix elements of processes (jHj)- (|15|) . excluding the Furry terms discussed below, are also 
implemented in the DISENT program package [T^j. Performing a numerical comparison 
with the latter, we find agreement. 

In Ref. ^H] j the squared matrix elements, which may be visualized as cut diagrams, are 
classified with respect to colour factors. One specific class, called F terms, contains all cut 
diagrams with two fermion loops, which are both coupled to three vector bosons, namely 
to one photon and two gluons. This class constitutes a gauge-parameter-independent sub- 
set of the NLO correction. The cut diagram of one specific member of this class is shown 
in Fig. El where the on-shell quarks are indicated by numbers. As was noticed in Ref. flB] , 
by Furry's theorem, each cut diagram within this class exactly cancels against one coun- 
terpart in which one fermion-number flow is reversed if the on-shell quarks associated with 
the loop whose fermion-number flow is reversed are not tagged in the experiment, i.e. if 
the three-momenta of these quarks are integrated over. This argument is also true in the 
case where only one fermion charge is identified, for instance in single-hadron production 
by e + e~ annihilation, since there is still a counterpart diagram where the other fermion- 
number flow is reversed. In our case, there are two tagged partons, one coming from the 
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proton and one fragmenting into the hadron h. Suppose the two tagged partons are the 
quarks 1 and 2 in the cut diagram of Fig. This situation can occur for processes (J12|) 
and (jl3)l . which involve only one quark flavour, and for process (I15J) . which involves two 
different quark flavours. Then, the Furry cancellation is impeded because there is no 
counterpart diagram. Thus, we are not allowed to omit this class of cut diagrams in our 
calculation. The corresponding squared matrix elements are listed in Appendix iBl 




Figure 2: Cut-diagram involving two fermion traces each coupled to three gauge bosons, 
as illustrated in Ref. [T%] . The cut proceeds along the numbered ticks representing the 
four on-shell quarks. If the charges of both quark loops are tagged, there is no Furry 
cancellation with an analogous diagram with one fermion-number flow flipped. 

The differential cross sections of processes (j8|)- ()15)) have to be integrated over the three- 
momenta of the second and third final-state partons keeping the three-momentum of the 
first one fixed. Performing the phase-space integrations, we encounter IR singularities of 
the soft and/or collinear types, which, for consistency with the virtual correction, must 
be extracted using dimensional regularization. It is convenient to do this by means of 
the dipole subtraction formalism jllj . The general idea of this formalism is to subtract 
from the contribution to the real correction due to a given partonic subprocess some 
artificial counterterm which has the same point-wise IR-singular behaviour in D space- 
time dimensions as the considered part of the real correction itself. Thus, the limit 
e — > can be performed, and the phase space integration can be evaluated numerically 
in four dimensions. The artificial counterterm is constructed in such a way that it can 
be integrated over the one-parton subspace analytically leading to poles in e. Adding 
the terms thus constructed to the virtual correction, the IR singularities of the latter 
are cancelled analytically. In the present case, where the three-momenta of two tagged 
partons need to be kept fixed, additional, more complicated artificial counterterms appear 
than in situations where only one parton is tagged, such as inclusive jet production. A 
technical advantage of the dipole subtraction method compared to the phase space slicing 
method is that all IR singularities cancel before any numerical integration is performed. 
Furthermore, there is no need to introduce a slicing parameter to separate soft and/or 
collinear phase space regions from the remaining hard region, which needs to be tuned 
in order to obtain a numerically stable result. For the factorization of the collinear 
singularities associated with the tagged partons, we choose the MS scheme. In turn, we 
have to employ PDFs and FFs which are defined in the same scheme. 



7 



Finally, we end up with two contributions, the real correction with the artificial coun- 
terterms subtracted and the virtual correction with the integrated artificial counterterms 
included, which are both finite in the physical limit e — > and can be integrated over 
their three- and two-particle phase spaces, respectively, in three spacial dimensions. These 
integrations are performed numerically using a custom-made C ++ routine. On the other 
hand, all algebraic calculations are executed with help of the symbolic-manipulation pack- 
age FORM Pl|. 

3 Numerical results 

We are now in a position to present our numerical results for the cross section of single- 
hadron inclusive production in ep DIS. We start by specifying our input. We work in the 
MS renormalization and factorization scheme with rif = 5 massless quark flavours. At 
NLO (LO), we employ set CTEQ6M (CTEQ6L1) of proton PDFs by the Coordinated 
Theoretical- Experimental Project on QCD (CTEQ) 2(1, the NLO (LO) set of FFs for 
light charged hadrons (7r ± , K^, and p/p) by Kniehl, Kramer, and Potter (KKP) \'21\ . and 
the two-loop (one-loop) formula for cts l! \p r ) with asymptotic scale parameter Aq^ d = 

226 MeV (165 MeV) [20J. This value is compatible with the result A^ D = (213±80) MeV 
((88 ±41) MeV) determined in Ref. [221- We approximate the ir° FFs as 

D?{x,ii f ) = ±D?{x,pf), (16) 

where D 1 ^ refers to the sum of the n + and n~ mesons, which is supported by LEP1 data 
of hadronic Z°-boson decays j2H|. Furthermore, we assume the charged hadrons to be 
exhausted by the charged pions, charged kaons, protons, and antiprotons, viz 

Djtfanf) = Df(x,p f )+Df(x,p f )+D^(x,p f ). (17) 

For simplicity, we identify the renormalization scale p r and the initial- and final-state 
factorization scales, pi and pf, respectively, and relate them to the characteristic dimen- 
sionful variables Q 2 and p* T by setting p? r — p 2 — p 2 — t;[Q 2 + (p^,) 2 ]/2, where £ is a 
dimensionless parameter of order unity introduced to estimate the theoretical uncertainty 
due to unphysical-scale variations. As usual, we consider variations of £ between 1/2 and 
2 about the default value 1. 

We now compare our theoretical predictions with HERA data on tt° mesons in the 
forward region from the HI Collaboration |3f4j and on charged hadrons in the current- 
jet region from the ZEUS Collaboration [3]. We start by discussing the HI data |3l4j . 
which were taken in DIS of positrons with energy E e = 27.6 GeV on protons with energy 
E p = 820 GeV in the laboratory frame, so that v^S* = 2-J ' E e E p = 301 GeV, during the 
running periods 1996 and 1996/1997, and correspond to integrated luminosities of 5.8 and 
21.2 pb _1 , respectively. In Refs. |3|4j . the ir° mesons were described by their transverse 
momentum p* T in the 7*p cm. frame and by their angle 9 with respect to the proton 
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flight direction, their pseudorapidity rj = — ln[tan(#/2)], and their energy E — x^Ep in 
the laboratory frame. They were detected within the acceptance cuts p^ > 2.5 GeV 
or 3.5 GeV, 5° < 9 < 25°, and xe > 0.01. The DIS phase space was restricted to 
the kinematic regime defined by 0.1 < y < 0.6 and 2 < Q 2 < 70 GeV 2 . The cross 
section was measured differentially in p^ |3|4| . 77 xe and xb |3|4| for various 
Q 2 intervals, differentially in for various xb intervals jl], and differentially in Q 2 [3]. 
The differential cross sections da n °/dp^, da'* /dr), da^/dxE, da n °/dxB, and da n °/dQ 2 
presented in Refs. [3] (open circles) and [3] (solid circles) are compared with our LO 
(dashed histograms) and NLO (solid histograms) predictions in Figs. 0H3 respectively. 
In Figs. El EJ Efa), andEfa), the upper three frames refer to the Q 2 intervals 2 < Q 2 < 
4.5 GeV 2 , 4.5 < Q 2 < 15 GeV 2 , and 15 < Q 2 < 70 GeV 2 . In Fig. EJb), the upper 
three frames refer to the x B intervals 0.000042 < x B < 0.0002, 0.0002 < x B < 0.001, 
and 0.001 < xb < 0.0063. In Fig. Efb), the upper three frames refer to the Q 2 intervals 
2 < Q 2 < 8 GeV 2 , 8 < Q 2 < 20 GeV 2 , and 20 < Q 2 < 70 GeV 2 . In all figures, the 
minimum-p^ cut is p^ > 2.5 GeV, expect for Fig. El^b), where it is p^ > 3.5 GeV. In 
Figs. 0HZI the shaded bands indicate the theoretical uncertainties of the NLO predictions 
due to the £ variation described above. The K factors, defined as the NLO to LO ratios 
of our default predictions, are shown in the downmost frames of Figs. 0HZI 

We observe from Figs.OHH that the HI data generally agree with our NLO predictions 
within errors, while they significantly overshoot our default LO predictions. Indeed, the K 
factors always exceed unity and even reach one order of magnitude at low values of p^, Q 2 , 
or xb- Not only do the LO predictions disagree with the HI data in their normalizations, 
but they also exhibit deviating shapes. On the other hand, under the effect of asymptotic 
freedom, the K factors approach unity for increasing values of fi r , i.e. for increasing values 
of p* T and/or Q 2 . 

There is an obvious explanation for the sizeable K factors at low values of \x r in terms 
of the different kinematic constraints at LO and NLO. The LO processes ©-(JTJ) are 
2^2, and their cross sections are sensitive to collinear singularities only as p* T — > 0. By 
contrast, processes dHJ (HSJ) contributing to the real NLO correction are 2 — > 3, so that 
collinear configurations can also arise for finite values of p* T . After mass factorization of 
the corresponding collinear singularities, the finite remainders can be sizeable, leading to 
large NLO corrections. A similar line of reasoning was presented in Ref. [12] . 

Unfortunately, the theoretical uncertainties in our NLO predictions due to £ varia- 
tion are rather sizeable, especially at low values of p* T , Q 2 , or xb, where the K factors 
themselves are abnormally large. This is partly related to the opening of new partonic 
production channels at NLO, which are still absent at LO, namely those of Eqs. ([lip. 
(JT3J), and ([To]) . Obviously, a reduction in £ dependence can only be expected to happen 
at next-to- next-to-leading order (NNLO), which is beyond the scope of this work. 

Besides the freedom in the choice of the renormalization and factorization scales, there 
are other sources of theoretical uncertainty, including the variations of the PDF and FF 
sets. However, in view of the considerable spread in cross section induced by the moderate 
£ variations described above, we conclude that the residual sources of theoretical uncer- 
tainty are of minor importance. Furthermore, we must bear in mind that the factorization 
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theorem itself is only valid up to terms of 0(Aq CD /(p^,) 2 ), which may become large in 
the \ow-p T range. 

We now turn to the ZEUS data on charged hadrons jHj, which were produced in DIS 
of electrons with energy E e = 26.7 GeV on protons with energy E p = 820 GeV in the 
laboratory frame, giving VS = 296 GeV, during the 1993 running period and correspond 
to an integrated luminosity of 0.55 pb™ 1 . They refer to the DIS phase space defined by 
10 < Q 2 < 160 GeV 2 and 75 < W < 175 GeV, where W is the 7*p invariant mass, with 
W 2 = (P + q) 2 = (1 — x B )yS, and come as multiplicities differential in p* T or Feynman's 
x variable Xp — 2p* L /W, where p* L = p^sinh^* is the projection of the hadron three- 
momentum onto the flight direction of the virtual photon in the 7*p cm. frame, and 
normalized to the total number of DIS events. Unfortunately, the xp distribution of 
the multiplicity includes charged hadrons with p^ values down to zero, while our NLO 
analysis is only valid for finite values of p^, so that a comparison is impossible. However, 
a comparison is feasible for the p* T distribution (l/N evt )dN ha _d/dp T 13 > which includes 
charged hadrons with xp > 0.05. The differential cross section dcr h± jdp* T may be obtained 
using the conversion formula 



1 do h 1 dNi 



had 



o"tot S dp* T N evt dp* T 
where cx°{. s is the total cross section in the DIS regime specified above, 



<t = dQ 2 dW———. (19) 



dQ 2 dW 



At LO, we have [2H| 

= W^x B [l + (1 - y) 2 }F2(x B , Q 2 ), (20) 
where xb = Q 2 /{Q 2 + W 2 ) and y = (Q 2 + W 2 )/S. Using the parameterization 

j \ c 2 +c 3 ln(l+Q 2 /Qo) 



F2(x B , Q 2 ) = Cl ^-j , (21) 

where Q 2 = 0.4 GeV 2 , c x = 0.2030±0.0086, c 2 = 0.0727±0.0046, and c 3 = 0.0448±0.0012, 
obtained from a fit to ZEUS data, we thus find <r°{ s = (35.4 ± 2.1) nb assuming the 
errors on c\, C2, and C3 to be statistically independent. This nicely agrees with the result 
a tot S = 33.9 nb obtained in the parton model of QCD, where j2S| 

F^x B ,Q 2 )=x B J2el [FV(x B ,Q 2 )+Fl(x B ,Q 2 )\ , (22) 

i=l 

using set CTEQ6L1 j20j of proton PDFs with rif = 5. For consistency, we use the 
ZEUS result for a°f to convert the ZEUS data for (l/N evt )dN hiid /dp* T 0. The result for 
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do h jdp* T thus obtained (solid circles) is compared with our LO (dashed histogram) and 
NLO (solid histogram) predictions in Fig. |H1 (upper frame). As in Figs. 0H3 the shaded 
band indicates the theoretical uncertainty in the NLO prediction due to the £ variation 
described above, and the K factor is also shown (lower frame). Again, our NLO prediction 
leads to a better description of data than our LO one. Here, the K factor takes more 
moderate values than under HI kinematic conditions, being of order 1.5 or below. As 
explained above, our LO and NLO predictions break down in the limit p* T — > 0. This 
drawback can be fixed by the resummation of multiple parton radiation, as demonstrated 
in Ref. |2E] on the basis of the LO result. 

In Section we explained why the Furry terms do not vanish in our case, in contrast 
to inclusive jet production in DIS ^J|- It is interesting to investigate their importance 
quantitatively. To this end, we reconsider the differential cross sections da n /dxs for 
0.1 < y < 0.6, 4.5 < Q 2 < 15 GeV 2 , p* T > 2.5 GeV, 5° < 6 < 25°, and x E > 0.01 and 
da h± /dp* T for 10 < Q 2 < 160 GeV 2 , 75 < W < 175 GeV, and x F > 0.05, which we already 
studied in the second frame of Fig. E^a) and the first frame of Fig. |HJ respectively, and 
turn off the Furry terms in our default NLO prediction. The results are shown together 
with our default LO and NLO predictions and the HI JT] and ZEUS [3] data in Figs. Efa) 
and (b), respectively. We observe that the Furry terms are very important. In Fig. Efa), 
they account for roughly 20% of the NLO correction, while, in Fig. E^b), they practically 
exhaust the latter. 

We expect our fixed-order predictions to break down in three extreme kinematic 
regimes corresponding to the limits (i) Q 2 — > 0; (ii) 6 — > or, equivalently, rj — > oo 
or xp — > — 1; and (iii) xb — > 0. Case (i) corresponds to the photoproduction limit, in 
which the resolved-photon contribution gains importance, especially at small values of 
and/or 9. Case (ii) is related to the possibility that the observed hadron h originates from 
the proton remnant, so that the notion of fracture functions is invoked. Case (iii) is ex- 
pected to correspond to the realm of BFKL |U] dynamics, although it is unclear precisely 
where the onset of the latter is supposed to be located. Our analysis is puristic in the sense 
that resolved virtual photons, fracture functions, and BFKL dynamics are disregarded, so 
as to test their actual relevance in the confrontation of the QCD-improved parton model 
with the experimental situation of Refs. |3|4|5j . Let us now scrutinize these issues. Doing 
this, however, we have to bear in mind that the theoretical uncertainty in our NLO pre- 
dictions due to the arbitrariness in the choice of the unphysical scales is particularly large 
in these corners of phase space, so that any conclusions are likely to be premature prior 
to the advent of a full NNLO analysis. From Fig. we observe that our NLO prediction 
for da n ° /dQ 2 tends to undershoot the HI data |3] in the low-Q 2 range, so that there is 
indeed some room for a resolved-photon contribution. Similar conclusions were reached 
in Ref. [27J. On the other hand, we see from Fig. 0] that the HI data for da n /drj [3J 
significantly exceed our NLO prediction in the very forward region, i.e. in the rightmost 
T) bin, for low values of Q 2 . In fact, for 2 < Q 2 < 4.5 GeV 2 , the measured rj distribution 
exhibits a plateau in the upper rj range, whereas the NLO prediction is rapidly suppressed 
by the shrinkage of the available phase space for increasing value of %]. This plateau might 
be partly caused by tt° mesons originating from the remnant jet, which contaminate the 
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proper data sample. Such events cannot be described within our puristic NLO QCD 
framework. Finally, thanks to the support from the Furry terms, we find in Figs. Efa) 
and (b) satisfactory overall agreement between our NLO prediction for da'" jdxs and the 
HI data [I] down to the lowest xb values considered. A similar conclusion can be drawn 
from Fig. Efb) for da w jdxE in the low-x^ bin 4.2 x 10~ 5 < xb < 2 x 1CT 4 . This suggest 
that, in the case of light-hadron inclusive production in DIS at HERA, the influence of 
the BFKL dynamics is likely to be still feeble for xb ^4.2 x 10 -5 . 

4 Conclusion 

We analytically calculated the cross section for the inclusive electroproduction of single 
hadrons with finite transverse momenta via virtual-photon exchange at NLO in the QCD- 
improved parton model, with rif massless quark flavours, on the basis of the collinear- 
factorization theorem. We worked in the MS renormalization and factorization scheme 
and handled the IR singularities using the dipole subtraction formalism [H]. As for the 
virtual correction, we reproduced the result of Ref. ^7j. As for the real correction, we es- 
tablished agreement with Ref. JSI, up to the Furry terms, which vanish upon phase-space 
integration in the case of single-jet inclusive electroproduction considered in Ref. [THJ, but 
yield a finite contribution in the case under consideration here. 

Using nonperturbative FFs recently extracted from data of e + e~ annihilation |21j . we 
provided theoretical predictions for the production of tc° mesons in the forward region 
and of charged hadrons in the current-jet region, and compared them in all possible ways 
with HI |3l4j and ZEUS [Sj data, respectively. Specifically, we considered cross section 
distributions in p^, r], xe, xb, and Q 2 . 

We found that our LO predictions always significantly fell short of the HERA data 
and often exhibited deviating shapes. However, the situation dramatically improved as we 
proceeded to NLO, where our default predictions, endowed with theoretical uncertainties 
estimated by moderate unphysical-scale variations, led to a satisfactory description of the 
HERA data in the preponderant part of the accessed phase space. In other words, we 
encountered K factors much in excess of unity, except towards the regime of asymptotic 
freedom characterized by large values of p* T and/or Q 2 . This was unavoidably accom- 
panied by considerable theoretical uncertainties. Both features suggest that a reliable 
interpretation of the HERA data |3l4l5j within the QCD-improved parton model ulti- 
mately necessitates a full NNLO analysis, which is presently out of reach, however. For 
the time being, we conclude that the successful comparison of the HERA data with our 
NLO predictions provides a useful test of the universality and the scaling violations of 
the FFs, which are guaranteed by the factorization theorem and are ruled by the DGLAP 
evolution equations, respectively. 

Significant deviations between the HERA data and our NLO predictions only occurred 
in certain corners of phase space, namely in the photoproduction limit Q 2 — > 0, where 
resolved virtual photons are expected to contribute, and in the limit rj — > oo, where 
fracture functions are supposed to enter the stage. Both refinements were not included 
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in our analysis. Interestingly, distinctive deviations could not be observed towards the 
lowest xb values probed, which indicates that the realm of BFKL j5] dynamics has not 
actually been accessed yet. 



Note added 



After finalizing this manuscript, a paper has appeared which also reports on a NLO 
analysis of the inclusive electroproduction of single hadrons with finite transverse momenta 
[2H] reaching conclusions similar to ours. 
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A LO results 

In this appendix, we list the LO expressions for Hj? and if£ b in Eq. (jlj) pertaining to 
processes ©-0, with ab = qq,qg,gq, respectively. We have 

H«? = lGnasCpe f + {1 ' X ~ ^ 



(1 -*)(!-*) 



Hf = %na s C F e 2 q Q 2 



H q T 9 = lQna s C F e 2 q l + { - X ^ 

Hl 9 = 87ra s C F elQ 2l ~ Z 
„ (/ 16na s N c C F e 2 g 1 - 2x(l - x) - 2z{l - z) 



•> (1 - x)z 
i 

19 — q^„, n_„2r>2- 



TT 

T ~ N£ — 1 z(l - z) 

i^flfflhf, (23) 

L N*-l x ' V ; 

where A^ c = 3 is the number of quark colours and C F = (N 2 — 1)/(2N C ) = 4/3 is the 
eigenvalue of the Casimir operator in the fundamental representation of the QCD gauge 
group SU(iV c ). 
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B Real correction: Furry terms 

In this appendix, we list the NLO expressions for Hf and Hf in Eq. flU) that originate 
from hindered Furry cancellations in the squared matrix elements of processes (II 2 j) and 
(|13|). with ab = qq, and of process (|T3|) . with ab = qq'. We denote the four- momenta of the 
second and third final-state quarks by p b and p c , respectively, and introduce the invariants 
Sij = Pi'Pj, where i, j = a, b, c, d with i ^ j. For given q, p a and pb, we need to integrate 
over p c , while pd is fixed through four-momentum conservation to be pd = q + p a — Pb — Pc- 
We work in the coordinate frame defined in Fig. ^ and parameterize p c as 



x c / Q 2 

Pr — —\ —, cos ol sin j3, sin a sin f3, cos/3), (24) 

2 y x(l — x) 

where a and /3 are the azimuthal and polar angles, respectively. Then, we have 

Ht,l > = F^a^bQ 2 / dx c j da j dcos(3h^, (25) 



where 



h F\qq _ _ L_^J s 2 _ S 2 ab S cd ~ S ac [s ad Sbc + S bd (s b d + S cd )] 

SSabScd (Pc- q) 2 I 

+ s ab [s ad s bc - 2s bc s bd + s ac {s bd - s cd ) - s bd Scd + s 2 cd ] 
- s a d{2s bd s c d + s bc (3s bd + Scd)] 



- (Pa <"> ~Pb) 
+ (Pb <-> Pc, Pa 

- (p C <"> Pd) 
+ (Pb Pc) 

- (Pa -> -Pc, Pb -> "Pa, Pc -> Pb) 
+ (Pa "Pd) 

- (Pb -> Pc, Pc -> Pd, Pd -> Pb), (26) 

^r, i 7 \ol v7 ^ 'ac^ bd($ a c ^ab^ad ^ac^ad) 

SabSacSbdScd I (Pb ~ q) (Pd ~ q) 2 



x ( s ad s bc 4" s ac s bd s ab s cd) 
+ 



^ab^cd^ab ^ab^ad $ ac$ ad) ab$ cd ~ SadSbc ~ °ac Sbd) 



(Pc - <?) 2 L(Pd - g) 

7 vjfSacSbd Sad^bc ~\~ S a bS c d){s ac S a dSbd{Sac ^ad) "I - ^ab^cdi^ac Sad 

(Pb - q) 2 

2 „ „2 



+ Sab(s ac S bd ~ S^Scd)) 



(27) 
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h: 




S 2 ac s bd ~ S 2 ab S cd - S ac [s ad S bc + S bd (s bd + S cd )] 



+ s ab [s ad s bc - 2s bc s bd + s ac (s bd - s cd ) - s bd s cd + s 2 cd ] 

- s ad [2s bd s cd + s bc (3s bd + Scd)] | 

- (Pa <-> "Pb) 

+ (Pb <"> Pc, Pa <-> "Pd) 

- (Pc <"> Pd) 



(28) 





2 [^a&^ac "I - S a d(s ac S a d)](s ac Sfcd S ad S bc -\- S ab S cd ) 



(29) 
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Figure 3: Differential cross section da 77 " /dp^ (in pb/GeV) of e + p — *> e + 7r° + X in DIS with 
0.1 < y < 0.6 and 2 < Q 2 < 4.5 GeV 2 (first frame), 4.5 < Q 2 < 15 GeV 2 (second frame), 
or 15 < Q 2 < 70 GeV 2 (third frame) at HERA with E e = 27.6 GeV and E p = 820 GeV 
for 7T° mesons with 5° < 9 < 25° and xe > 0.01. HI data from Refs. [Sj (open circles) and 
[3] (solid circles) are compared with our default LO (dashed histograms) and NLO (solid 
histograms) predictions including theoretical uncertainties due to £ variation (shaded 
bands). The K factors (fourth frame) are also shown. 
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Figure 4: Same as in Fig. El but for da n /drj (in pb) with p* T > 2.5 GeV. 



18 



10 

1 

0.1 



r-i 0.01 

■§ 10 



X 



0.1 



0.01 
10 



1 

0.1 

0.01 
10 

o 



2 <Q"<4.5GeV 



5 — 



.4A<QI<J_5_GeV_ 
15<Q 2 <70GeV 2 



0.02 



~> — 2 1 r 

2 < Q < 4.5 GeV 



• HI data '04 
-- LO 
NLO 

i 



4.5 < Q 2 < 15 GeV 2 



L I 



15 < Q 2 < 70 GeV 2 



0.04 



0.06 



Figure 5: (a) Same as in Fig. but for da n /dx E (in nb) with p* T > 2.5 GeV. 
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Figured (b) Same as in Fig.Efa), but for 2 < Q 2 < 70 GeV 2 and 0.00042 <x B < 0.0002 
(first frame), 0.0002 < x B < 0.001 (second frame), or 0.001 < x B < 0.0063 (third frame). 
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Figure 6: (a) Same as in Fig. El but for da 7 " /dx B (in nb) with p* T > 2.5 GeV. 
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Figure EJ (b) Same as in Fig. EJa), but for p* T > 3.5 GeV and 2 < Q 2 < 8 GeV 2 (first 
frame), 8 < Q 2 < 20 GeV 2 (second frame), or 20 < Q 2 < 70 GeV 2 (third frame). 
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Figure 7: Same as in Fig. El but for da" /dQ 2 (in pb/GeV 2 ) with p* T > 2.5 GeV. 
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Figure 8: Differential cross section da h /dp^ (in pb/GeV) of e p — > e + X in DIS 
with 10 < Q 2 < 160 GeV 2 and 75 < W < 175 GeV at HERA with E e = 26.7 GeV and 
E p = 820 GeV for charged hadrons with xp > 0.05. ZEUS data [5] (solid circles) are 
compared with our default LO (dashed histograms) and NLO (solid histograms) predic- 
tions including theoretical uncertainties due to £ variation (shaded bands). The K factor 
(lower frame) is also shown. 
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Figure 9: Same as in (a) the second frame of Fig. IHfa) and (b) the first frame of Fig. |Hl 
but also including our default NLO predictions with the Furry terms turned off (dotted 
histograms). For clarity, the theoretical uncertainties due to £ variation are omitted. 
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